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*n ABSTRACT 

o 

Medium-resolution spectra from 3650 A to 10,000 A are presented for 96 giant H II 
regions distributed in 20 spiral galaxies. In order to interpret the data, we have cal- 
culated two separate grids of photoionization models, adopting single-star atmospheres 
(Kurucz) and star clusters synthesized with different Initial Mass Functions (IMFs) as 
ionizing sources. Additional models were computed with more recent non-LTE stellar 
atmospheres, in order to check the effects of different stellar ionizing fluxes. We use 



Of 
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the radiation softness parameter rj/ = ([OII]/[OIII])/([SII]/[SIII]) of Vilchez and Pagel 
to test for a metallicity dependence of the effective temperatures of the ionizing stars. 
Our results are consistent with a significant decrease in mean stellar temperatures of 
the ionizing stars with increasing metallicity. The magnitude of the effect, combined 
with the behavior of the He I A5876/H/? ratio, suggest a smaller upper mass limit for 
star formation at abundances higher than solar, even when considering the effects of 
metallicity on stellar evolution and atmospheric line blanketing. However, the exact 
magnitudes of the stellar temperature and IMF variations are dependent on the choice 
of stellar atmosphere and evolution models used, as well as on uncertainties in the 
nebular abundance scale at high metallicities. 

Our results also constrain the systematic behavior of the ionization parameter and 
the N/O ratio in extragalactic H II regions. The observed spectral sequences are in- 
consistent with current stellar evolution models which predict a luminous, hot W-R 
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stellar population in evolved H II regions older than 2-3 Myr. This suggests either 
that the hardness of the emitted Lyman continuum spectrum has been overestimated 
in the models, or that some mechanism disrupts the H II regions before the W-R phases 
become important. 

Subject headings: galaxies: ISM — galaxies: spiral 

1. INTRODUCTION 

Giant extragalactic H II regions are ionized by the massive stars of their young embedded clus- 
ters, and as such provide powerful tools for measuring the chemical composition of the interstellar 
medium and for understanding the massive stellar populations in external galaxies (see reviews 
by Shields 1990, Kennicutt 1991, Stasihska 1996). Evidence for abundance gradients across the 
disks of spiral galaxies was presented by Searle (1971), following pioneering work by Aller (1942) 
on the radial trends of emission-line ratios. The excitation of H II regions, measured by the ra- 
tios [OIII]AA4959,5007/H/3 and [OIII]A5007/[OII]A3727, was found to increase with galactocentric 
distance, and was interpreted as an effect of the shift from infrared fine-structure transitions to op- 
tical [OIII] emission as the principal cooling mechanism at the lower metallicities typical of larger 
radii. Later work confirmed the existence of composition gradients by direct observation of the 
[OIII]A4363 line at low abundances (e.g. Smith 1975). To explain the observed gradients in ex- 
citation and in the equivalent width of the H/3 emission [EW(H/?)] an additional gradient of the 
color temperature of the ionizing radiation field, T*, was proposed (Shields 1974, Shields & Tinsley 
1976, Shields & Searle 1978), with the ionizing radiation becoming harder with increasing radius. 
The work of Shields & Tinsley (1976) suggested further that this change in mean stellar temper- 
ature was produced by a metallicity-dependent upper mass limit for star formation, M u oc Z -1 / 2 . 
The same authors also pointed out that a change in the temperature of the hottest stars is to be 
expected, even in the absence of a metallicity effect on M u , because of opacity effects that would 
lower the effective temperatures of metal-rich stars. 

Subsequent efforts were concentrated on quantifying this general picture, with the calculation 
of extensive grids of single-star photoionization models. Most authors confirmed the increase of 
T* with decreasing abundance (Stasihska 1980, Campbell, Terlevich & Melnick 1986, Vflchez & 
Pagel 1988), together with a possible increase in the ionization parameter (Campbell 1988). On 
the other hand Evans & Dopita (1985) argued that the available data on the trends in the emission- 
line ratios were consistent with a single ionizing temperature (~ 41,500 K) and a varying nebular 
geometry and ionization parameter. As the controversy is not completely settled, the presence of 
a gradient still needs to be clearly established and quantified, over the full range of abundances 
observed in spirals. 

The fundamental question remains whether these apparent changes in T* are due to systematic 
variations in the Initial Mass Function (IMF) with metallicity and/or galaxy type. Abundance 
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effects on either the upper mass limit (Shields & Tinsley 1976) or the slope of the IMF (Terlevich 
1985) have been proposed. However, several other mechanisms have been suggested, including 
increased stellar line blanketing (Balick & Sneden 1976) or internal dust absorption (Sarazin 1976) 
with metallicity. Modern stellar evolution models (e.g., Maeder 1990) predict a softening of the 
stellar ionizing spectra with increasing abundance. The consequence for H II region models is a 
natural decrease of T* , making the varying IMF hypothesis less compelling, as shown by McGaugh 
(1991). Moreover, the introduction of cluster photoionization models, in which the single star 
is replaced by an evolving stellar cluster as the ionizing source, has produced results which are 
consistent with a universal IMF (Garcia- Vargas, Bressan & Diaz 1995, Stasihska & Leitherer 1996). 

In this paper we re-examine the evidence for T* and IMF gradients, using new spectra of 96 
H II regions in 20 galaxies over the 3650-10,000 A range. This work extends previous large surveys 
of H II region spectra in spiral disks, including those of McCall, Rybski & Shields (1985), Vila- 
Costas & Edmunds (1992), Oey & Kennicutt (1993) and Zaritsky, Kennicutt & Huchra (1994). 
The main addition of the present work is the inclusion of data for the near-IR [SIII]AA9069,9532 
lines, which allows measurement of the radiation softness parameter rj/, directly related to the 
ionic fraction ratio rj = (0 + /0 ++ )/(S + /S ++ ), introduced by Vflchez and Pagel (1988) as a %g 
sensitive parameter. This paper extends to a large H II region sample the rjf method, which so 
far has been applied to studies of a few H II regions in individual galaxies. When combined with 
the empirical abundance indicator i?23= ([Oil] + [OIII])/H/3 (Pagel etal. 1979) and an extensive 
set of nebular models we can simultaneously constrain the metallicity, the effective temperature 
and the ionization parameter (Diaz etal. 1987, Vflchez etal. 1988, Diaz etal. 1991). A companion 
paper presents detailed modeling of 40 H II regions in one galaxy (M101), with emphasis on the 
T* gradient in a homogeneous sample of H II regions (Garnett & Kennicutt 1998). In the present 
work we examine a larger galaxy sample, spanning a wide range of abundance, with two main 
goals: to study the relationship bewteen T* and the gaseous abundance, and to derive constraints 
on IMF variations with metallicity. For this purpose, photoionization models including as ionizing 
source either an evolving star cluster or a single stellar atmosphere are computed. The new H II 
region sample is presented in § 2, and the nebular model grid is described in § 3. In the following 
sections we discuss our main results: systematics of H II region spectra and physical conditions 
(§ 4), evidence for a T* gradient (§ 5), and consequences for the upper IMF (§ 6). We summarize 
the main results in § 7. As a quick reference a list of abbreviations and symbols used is provided 
in Table |. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Sample 

The galaxy sample (Table [2J) spans the Hubble sequence from type Sa to type Sm, and should 
therefore provide good indication of trends of star forming region properties with morphological 
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type. The H II regions cover a wide range in excitation, and their chemical abundance varies 
from Z ~ 0.2 Z© to above solar (the lower limit is set by abundance determinations found in the 
literature for a few of the most metal-poor objects). Nearly all of the observations discussed here 
are new to this work, and were obtained during several runs in 1988-1992 on the Multiple Mirror 
Telescope (MMT) and the Steward Observatory 2.3m Telescope; the few exceptions are indicated 
in Table @. The main observational material (redwards of H/3) consists of spectra obtained with 
the MMT Red Channel Spectrograph with an echellette grating. Most of the blue spectra (which 
include the [Oil] A3727 line) were obtained at the MMT as part of the survey by Oey & Kennicutt 
(1993). Blue data for a few individual objects were drawn from the literature (McCall etal. 1985, 
Zaritsky etal. 1994, and Dufour etal. 1980). 



2.2. MMT Echellette Spectra 

The Red Channel Spectrograph at the MMT was used in cross-dispersed mode to obtain 
spectral coverage from 4700 A to 10,000 A. This allowed the inclusion of all the main spectral 
features from H/3 to [SIII] AA9069, 9532. In this observing mode a 2" x 20" slit was cross-dispersed 
into five orders, providing extended wavelength coverage at a resolution of 7-15 A. The detector 
used was a Texas Instruments thinned 800x800 device which, when binned 2x2, provided a spatial 
scale of / .'625 pixel -1 . Exposure times ranged from 5 to 40 min. 

Data reduction included a few procedures particular to the echellette observations (see also 
Skillman & Kennicutt 1993 and Kennicutt & Garnett 1996). Standard star frames taken at low 
airmass were used to trace the aperture for each of the five orders and to remove the curvature in- 
troduced by the cross-disperser. Wavelength and flux calibrations were accomplished with HeNeAr 
lamp frames and observations of standard stars from Oke & Gunn (1983) and Massey etal. (1988), 
respectively. The calibrations were checked for internal consistency using lines in the overlap re- 
gions between adjacent orders. The residuals from the standard star fits indicate a typical accuracy 
of ±0.05 mag rms or better in all orders except the reddest one (A > 8000 A). Here the water 
absorption features in the near-IR introduce a larger uncertainty in the flux calibration. The ex- 
tinction was then modeled as a function of wavelength by tracing the absorption features observed 
in the standard star spectra. This procedure satisfactorily removes most of the unwanted telluric 
features around the [SIII] AA9069, 9532 lines, as indicated by the mean observed A9532/A9069 ratio 
= 2.47 ± 0.05 (the theoretical value is 2.44). Nevertheless we attach a conservative uncertainty of 
±20% to the [SIII] fluxes (see Kennicutt & Garnett 1996). One-dimensional spectra were extracted 
with variable apertures, typically eight arcseconds, but dependent on the seeing and the size of 
the objects, centered on the peak of the Ha emission line. Care was taken to match as closely as 
possible the regions sampled in the blue spectra. 



- 5 - 



2.3. Blue Spectra 

Coverage of the 3650-5100 A spectral range was obtained with two different instruments and 
setups. The majority of our targets in early-type galaxies were observed with the Red Channel CCD 
spectrograph on the MMT and a slit size of 2" x 180". These observations have been described in 
detail by Oey & Kennicutt (1993) in their project on abundances in early-type galaxies. However, 
new one-dimensional extractions were made, to match those of the red echellette spectra. 

Additional objects (marked 90" in Table |2|) were observed with the B&C spectrograph at the 
Steward Observatory 2.3m telescope. Observations were made with a 600 groove mm -1 grating 
blazed at 3568 A, providing 9 A resolution with a 4"5 x 180" slit. Standard reduction procedures 
for these long-slit data were followed. The H II regions in three galaxies (M31, M33 and NGC 2403) 
were observed in the red using a 300 groove mm -1 grating blazed at 6690 A, and providing 20 A 
resolution with the 4"5 slit. These data provided a check on the corresponding MMT echellette 
spectra. Blue data for NGC 628, NGC 2403, NGC 3521 and NGC 5236 were taken from the 
literature, with references given in Table |2[ 

2.4. Analysis of Spectra 

The emission line fluxes for the atomic transitions of interest were measured with the IRAFQ 
task SPLOT, by integrating over the line profiles. We assumed a stellar underlying continuum 
Balmer absorption equivalent width of 2 A, the mean value found by McCall etal. (1985) in their 
H II region sample, and calculated the reddening with the Ha/H/3 and H/3/H7 ratios, assuming Case 
B recombination ratios for an electron temperature of 10,000 K and density of 100 cm~ 3 (Osterbrock 
1989). Table lists the measured fluxes, corrected for reddening, for the most important emission 
lines, as well as offsets in arcseconds from the galaxy nucleus. Uncertainties in line fluxes were 
derived following the procedure described in Kennicutt & Garnett (1996), who analyzed H II 
region spectra in M101 taken with the same setup as the observations discussed in this work. The 
quoted errors were calculated by adding in quadrature the contributions from the statistical noise, 
the calibration (±3%), the flat fielding (±1%), the relative scaling of the echellette orders (±4%) 
and the uncertainty in the reddening corrections. A 20% uncertainty was assigned to the combined 
[SIII] AA9069, 9532 flux, as described earlier. 

We performed a few tests in order to check the quality of the data. As mentioned above the 
observed values of the [SIII] A9532/ [SIII] A9069 ratio are consistent with the theoretical expectations. 
The density-sensitive [SII]A6716/[SII]A6731 ratios fall near or below the low-density limit (= 1.42), 
with an average value of 1.37 ± 0.02, therefore justifying the assumption of low density (10-100 
cm -3 ) made for the model nebulae (c/. § ||). Finally, in Figure [l] our data (filled circles) are 



3 IRAF is distributed by the National Optical Astronomical Observatories, which are operated by AURA, Inc., 
under contract to the NSF. 
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compared to the data of McCall etal. (1985, open circles) for two excitation sequences, [OII]/H/3 vs 
[OIII]/H/3 and [OIII]/H/3 vs [Nil] /Ha. As can be seen, the sequences for the two datasets closely 
follow each other. 



3. PHOTOIONIZATION MODELS 

In order to interpret the emission-line data for our H II regions, a set of nebular models was 
calculated with the CLOUDY program (Ferland etal. 1998; Ferland 1996), version 90.01. Details 
about the models and more extensive results can be found in Bresolin (1997). 

Given the spectral energy distribution (SED) of the ionizing source, the ionizing photon lu- 
minosity, the chemical composition and the geometry of the gas, the code solves for the ionization 
and thermal equilibrium of the nebula. A relevant aspect of the version of the code used for these 
calculations is the inclusion of Opacity Project data (Seaton etal. 1994). For our models the im- 
portant change (approaching a factor of two) with respect to previous versions of the code regards 
the + photoionization cross section, which no longer shows a jump at the 2s - 2p edge. A further 
change in the + - ++ balance is caused by new collision strengths for some ++ lines. To assess 
the effect of the new atomic data on the final results, a set of single-star photoionization models 
(described in the next section) were calculated with an older version of CLOUDY {v. 84.12), and 
compared with the results of v. 90.01. The newer version predicts an average decrease of the [Oil] 
line intensity of approximately 20%, mostly counterbalanced by an increase in the [OIII] emission. 
The corresponding effect on the empirical abundance indicator R23 translates into an approximate 
0.025 dex increase in the (O/H) abundance scale, based on the calibration given by Zaritsky etal. 
(1994). 

For the remainder of this section the input parameters of the models will be described; the 
relevant model results will be illustrated in § |I] 

3.1. Input Parameters 

3.1.1. Spectral Energy Distributions 

As mentioned above, two different SEDs were considered as the ionizing source of the nebulae. 
The single-star models are simpler cases in principle, as they can be interpreted in terms of one 
parameter, T*, which is a function of the hardness of the output ionizing spectrum. On the other 
hand, cluster models contain information on the mass function and the temporal evolution of the 
ionizing clusters. 

For single-star models, the line-blanketed LTE Kurucz (1992) stellar atmospheres, which in- 
clude the effects of bound-free opacity from heavy elements, were adopted for consistency with 
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the cluster models. The logg = 5 atmospheres were chosen, in order to circumvent problems of 
non-convergence of the hydrostatic LTE atmospheres at low stellar surface gravity. The metallicity 
was varied in accordance with the gas abundance. Nebular models were calculated at Z = 0.25 Z , 
Z and 2 Z , and for % = 37,500 K, 40,000 K, 45,000 K, 50,000 K and 60,000 K. The sensitivity 
of the results to the choice of stellar atmosphere models will be discussed later. 

An alternative approach is to calculate SEDs for star clusters with different IMFs. This 
method was applied by McGaugh (1991) for zero-age clusters, and taking into account the temporal 
evolution of the clusters by Garcia- Vargas & Diaz (1994), Cerviho & Mas-Hesse (1994), Garcia- 
Vargas etal. (1995), and Leitherer & Heckman (1995). 

We adopted the SEDs calculated by Leitherer & Heckman (1995), as presented in Leitherer 
etal. (1996a). These SEDs were also calculated for different metallicities (those at 0.25 Z , Z and 
2 Z are relevant to the present work) and IMFs (a = 2.35 and M up = 100 M , a = 2.35 and 
M up = 30 M , a = 3.35 and M up = 100 M , where the IMF is given by dN/dm oc mr a between 
Mzou> = 1 M and M up ), at 1 Myr intervals from 1 to 25 Myr. The stellar models of Maeder 
(1990) and Maeder & Meynet (1988) were used for the synthesis, together with Kurucz's (1992) 
atmospheres, supplemented by the non-LTE models of Schmutz, Leitherer & Gruenwald (1992) for 
stars with strong winds. We used the results of both the instantaneous burst case, in which star 
formation occurs over a timescale short compared to the age of the model, and the continuous star 
formation case (star formation constant with time). We considered ages from 1 to 6 Myr, since at 
later times most clusters would lack the necessary ionizing luminosity to produce an H II region 
(Garcia- Vargas & Diaz 1994). For further details on the cluster SEDs, the reader is referred to the 
original papers. 



3.1.2. Nebular Gas Parameters 

The abundance of the gas has been made consistent with the metallicity of the ionizing stars. 
The solar elemental abundances are adopted as defined in Grevesse & Anders (1989). To take into 
account depletion of refractory elements onto dust grains, the elements Mg, Al, Ca, Fe, Ni and Na 
were depleted by a factor of 10, and Si by a factor of two, relative to the solar abundance (consistent 
with Garnett etal. 1995, Garnett & Kennicutt 1998). The He abundance was scaled according to 
Y = Y p + (AY/AZ) Z, where Y = 4y/(l + Ay), y = He/H by number, and Z is the metallicity 
(Z = 0.02). We have adopted Y p = 0.23 (Pagel etal. 1992), and AY/AZ=2.5. Table | summarizes 
the adopted solar composition, with the above-mentioned depletions taken into account. 

The model geometry is spherical, with a constant gas density n = 50 cm" 3 , and with a filling 
factor e adjusted so as to obtain models with three different ionization parameters, log U = —2, —3 
and —4. The dimensionless ionization parameter U is defined as 



(1) 
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where Q H o is the number of hydrogen-ionizing photons (E > 13.6 eV) emitted per second, R s is 
the Stromgren radius, n is the hydrogen density and c is the speed of light. U depends on the 
electron temperature, T e , through the Stromgren radius dependence on the case B recombination 
coefficient ob(T): 



with A oc T e . We had therefore to iterate a few times, from an initial guess for the filling factor, 
to obtain the desired log?7. A few models were calculated for n = 10 cm~ 3 and n = 150 cm~ 3 , to 
estimate gas density effects, which are quite strong at higher metallicities (Diaz etal. 1991). 

We fixed the total number of ionizing photons to Q^o = 10 51 s" 1 , the equivalent of ~90 07 V 
stars (Vacca 1994), and representative of bright extragalactic H II regions. The exact value of Quo 
is of no practical importance, since it acts as a scaling factor, and models with a given ionizing 
spectrum shape and differing Q H o , n and e, but with the same U are homologous, producing the 
same emission-line relative strengths (as long as the density is below the critical limit for collisional 
de-excitation) and the same nebular ionization structure. The corresponding initial cluster masses 
are large enough (> 1 x 10 4 M ) that our models remain unaffected by the stochastic effects 
described by Garcia- Vargas & Diaz (1994) and Cerviho & Mas-Hesse (1994). We also emphasize 
that CLOUDY computes nebular models in which a compact ionizing source lies at the center of 
a spherically symmetric gas distribution. We do not expect significant changes in the case of a 
more extended ionizing source (e.g. a loose OB association) as long as it remains unresolved by 
the spectrograph slit, since the emission-line properties are controlled by a quantity integrated over 
the volume of the nebula (the ionization parameter). 

All of our models are ionization-bounded, with the calculation stopped at a gas temperature of 
100 K (the optical opacity of the gas is negligible even at a much higher temperature). The models 
include a 30 km s _1 turbulent velocity term, consistent with observed linewidths of giant H II 
regions (Arsenault & Roy 1988). This additional velocity field increases the line widths (already 
broadened by thermal motions), affecting the energy budget of the H II regions. A cooling of the 
nebulae results through infrared fine-structure lines, which can become optically thin as a result of 
the line broadening. The effects of dust have not been taken into account, except for the depletion 
of refractory elements. 




(2) 



with a B (T) oc T" 1 , 



so that we can rewrite 



U = A(Q H one 2 ) 1/3 



(3) 



For the interested readers, the results of our models can be retrieved, in the form of tables and 
diagnostic diagrams, by anonymous ftp, by contacting the first author. 
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4. RESULTS 

A useful quantity that can be readily estimated from the stellar atmosphere SEDs is the ratio 
of He ionizing photons (A < 504 A) to H ionizing photons (A < 912 A), Que IQh°i f° r each 
metallicity and %g. For a particular atmosphere model, this ratio describes the hardness of the 
radiation field, and has been used to define an equivalent ionizing temperature, T e<? , as the %g of 
a star whose SED produces a given QHe /QH (Mas-Hesse & Kunth 1991, Garcia- Vargas & Diaz 
1994). A calibration of this ratio for the Kurucz (1992) log g = 5 atmospheres is given in Table 
Clearly the radiation field becomes harder as T^jj increases, and with lowering metallicity at a fixed 
temperature. The sensitivity to abundance is largest at high temperatures, becoming progressively 
smaller at lower The abundance effects can be quite important, since they can be equivalent 
to a change in T* of as much as 10,000 K between 2 Z and 0.1 Z . 

The Qne° /Qh° ratios for the cluster models are shown in Table ||. As pointed out by Garcfa- 
Vargas etal. (1995), the evolution of T eq is non-monotonic at higher metallicities, due to the 
hardening of the radiation after the predicted W-R phase sets in. This effect disrupts the T eq - 
metallicity relationship valid for ages smaller than about 3 Myr, for which T eq drops by 5000-6000 K 
(from 50,000 to 45,000 K at t = 1 Myr) between 0.1 Z and 2 Z . The presence of the W-R stars 
strongly affects the results 2-3 Myr after the burst of star formation. The effect is stronger with 
increasing metallicities, because of the larger number of W-R stars produced. 

A direct comparison of our model results with those of other authors is made difficult by the use 
of different input parameters and photoionization codes. However, a small set of solar composition 
models, with n = 10 cm~ 3 , was compared to the models by Stasihska & Leitherer (1996). The 
latter were calculated with the code PHOTO (Stasihska 1990) and the same cluster SEDs used in 
the present work, for different cluster masses. Due to the evolution of the massive stars, the total 
ionizing luminosity (and hence the ionization parameter) is a function of cluster age, while in our 
models we have kept it constant, in order to set the ionization parameter to the desired values. The 
Stasihska & Leitherer models for cluster masses M* = 10 3 and 10 6 M and ages between 1 and 6 
Myr (M up = 100 M ) have an ionization parameter in the range -1.6 to -2.9, and can therefore be 
compared to our log U = — 2 and -3 models. Figure || shows the comparison regarding [OIII]A5007 
vs H/3 and [OII]A3727 vs H/3. There is a good agreement, considering the differences in ionization 
parameter between the two sets of models, and the different photoionization codes used. 

4.1. Ionization Parameter 

The model results indicate that the ionization parameters adopted for the models well bracket 
our data, which are clustered around log U = —3. This is demonstrated in Figure ||, which plots 
[SII]AA6717,6731/[SIII]AA9069,9532 vs the abundance parameter R23 for our sample. Diaz etal. 
(1991) gave a calibration of U in terms of this ratio, which has been refined by Garnett Sz Kennicutt 
(1998) to include the small metallicity dependence. Similarly we have used a surface fit to our 
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models in order to derive the following relation: 



log U = (-3.35 ± 0.04) - (1.54 ± 0.05) log {[SII]/ [SIII]) - (0.08 ± 0.06) log(i? 23 ) 



(4) 



According to this calibration and the observed [SII] /[SIII] ratios the ionization parameter in 
our sample is in the range \ogU = —3 ± 0.6. The comparison of the models with the observations 
suggests also that for the most metal-poor objects (~ 0.2 Z©) the ionization parameter is about 4 
times larger than for H II regions around the solar metallicity. This finding is in agreement, at least 
qualitatively, with previous investigations of H II galaxies (Campbell 1988, Stasihska & Leitherer 



The [Nil] /[Oil] model sequences do not reproduce the observations at low abundances when 
a constant N/O ratio is adopted. The actual data show weaker [Nil] lines relative to [Oil] when 
compared to the theoretical values. Theory and observations could be reconciled by considering a 
metallicity-dependent N/O ratio. Evidence for such an effect in spirals has already been shown by 
several authors {e.g. McCall etal. 1985, Fierro, Torres-Peimbert & Peimbert 1986, Diaz etal. 1991, 
Garnett & Kennicutt 1998), but there is no general agreement on the functional relation between 
N/O and O/H, if it exists at all (different spirals having possibly different N/O gradients, Henry & 
Howard 1995; these authors, however, neglected primary nitrogen in deriving their gradients). In 
low-abundance irregular galaxies the N/O ratio is found to be independent of O/H (Garnett 1990), 
supporting the idea of a primary origin for nitrogen at low metallicity: nitrogen is produced out of 
C and O synthesized during the evolution of stars. An abundance-sensitive N/O ratio would imply 
the presence of a secondary component, that is N produced by C and O already present in stars 
at their birth. The situation is less clear for S/O, for which the existence of a O/H dependence is 
uncertain (for example, see contradicting results in Diaz etal. 1991 and Garnett 1989). To represent 
a secondary component superposed on a primary component for N in spiral galaxies, a few models 
were calculated assuming the following relations (Garnett & Kennicutt 1998): 



1996). 



4.2. The N/O abundance ratio 



log N/O 



1.5 + (log O/H + 3.7) 



QogO/H > -3.7) 



(5a) 



log N/O = -1.5 



(log O/H < -3.7) 



(5b) 



The results (Fig. |j) are in much better agreement with the observations, lending more credibil- 
ity to the presence of a secondary component for nitrogen in spiral galaxies. The models can now 
reproduce the steepening of the [Nil] /[Oil] vs R23 relation that occurs in the data at Log R23 > 0.5. 
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The main secondary effect of the abundance-dependent N/0 ratio is a shift of the 2 Z models to 
lower R23 values, due to the increased cooling of the nebulae caused by the larger abundance frac- 
tion of nitrogen. Models with a fixed, solar N/O ratio would therefore overestimate the metallicity 
for high-abundance objects. The models presented in this paper were calculated with a constant 
N/O ratio; however the effects of a variable ratio for all the diagnostics but those involving [Nil] 
are, for all our purposes, marginal and limited to Z > Z . 

An important point to make about the [Nil] /[Oil] sequence is that it refers to 20 spiral galaxies 
of different Hubble type. The tightness of the sequence implies that the relation between N/O and 
O/H is similar for all these galaxies, whatever the exact functional form might be. While for a 
few individual galaxies (e.g. M33, see Vflchez etal. 1988, or M81, see Garnett & Shields 1987) a 
N/O gradient seems not to be required, a 'universal' abundance-dependent N/O ratio seems to be 
necessary to explain the bulk of our observations. 

4.3. The H/3 Equivalent Width 

Bresolin & Kennicutt (1997) studied the distribution of the Ha equivalent width [EW(Ha)] of 
H II regions in ten galaxies of different Hubble type to test the idea that different upper mass limits 
of the IMF could be responsible for the changes in star forming region properties as a function of 
morphological type. The photometrically-determined EW(Ha) values showed no dependence on 
Hubble type. The distribution of the equivalent width of Ha and H/3 for our H II region sample 
confirms this result. Since the equivalent width of the Balmer emission lines is a measure of the 
ratio of the number of hot, massive stars (producing the emission line flux through ionization of 
the nebula) to the total number of stars (producing the continuum), our result points toward the 
invariance of the IMF across the morphological sequence. 

The data show a dependence of EW(H/3) on metallicity, with a large scatter, likely due to an 
age spread in the H II region sample. A similar behavior was already measured by Searle (1971) 
on a small sample of H II regions in M101. This result was interpreted by Shields & Tinsley (1976) 
as evidence for a gradient in the temperature of the hottest exciting stars, which in turn led them 
to hypothesize an abundance-dependent upper mass limit for star formation. The photoionization 
models, however, fail to reproduce the observed sequence of EW(H/3) vs R23. It is a well-known fact 
that the observed EW(H/3) values of extragalactic H II regions are much lower than the theoretical 
predictions. The usual interpretation is that an underlying population of stars older than the ones 
responsible for the ionization of the nebulae is contributing to the stellar continuum, thus diluting 
the ionizing radiation and reducing the equivalent width (McCall etal. 1985). Diaz etal. (1991) 
proposed the coexistence of stellar clusters of different ages within giant H II regions, as suggested 
by spatially resolved observations (Skillman 1985, Diaz etal. 1987). Internal extinction due to dust 
which would affect the nebular lines, but not the continuum, has been proposed as an alternate 
solution (Mayya & Prabhu 1996, Garcia- Vargas etal. 1997). The dependence of EW(H/3) on both 
cluster age and upper IMF, combined with the problem just described, make it difficult to derive a 
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leff-'Z relationship from EW(H/3). Therefore we will not attempt to use the observed EW(H/3) to 
infer a T* gradient. 

4.4. Cluster Ages 

A T*-metailicity relation is expected for young clusters (1-2 Myr after initial burst), as a result 
of metallicity-dependent effects on the stellar output flux. The onset of the W-R phase for massive 
stars which occurs after 2-3 Myr of the initial burst would disrupt this correlation for a mixed-age 
sample of clusters, as shown by Garcia- Vargas etal. (1995). The observation of a T* gradient in 
our sample of H II regions (§ 5) would therefore imply, if theory and models are correct, that an 
age limit exists. Interestingly, this would agree with the fact that the instantaneous burst cluster 
models are only able to reproduce the data for young cluster ages (t < 3 Myr). The theoretical 
predictions diverge considerably from the observations at later times. This does not exclude the 
likely presence of older H II regions in our sample, but the average of the H II region sequences 
is best reproduced by the 1 and 2 Myr models, when we make the assumption that the ionizing 
sources are single population clusters. As an example, we show the diagnostic diagrams which 
summarize the model results for [OII]/[OIII] vs R23 in Figure |H[ Here the data (solid dots) are 
compared with the models, calculated for three different values of the ionization parameter and 
for three different cluster IMFs, as explained in § ^. It can be seen that the hardening of the 
stellar radiation predicted after 2 Myr from the initial burst produces line ratios which cannot in 
general be reconciled with the observations (see also Fig. |l3| ). There remains the possibility that 
the nebular models are at fault, even though state-of-the-art codes and models have been adopted 
for both the gaseous and the ionizing source treatment. Alternatively some of the predictions of the 
stellar evolutionary models could be incorrect. In short, our results imply that either the current 
evolution models overpredict the hardness of the emitted Lyman continuum spectrum, or that some 
physical mechanism, such as disruption by stellar winds and/or supernova explosions, prevents us 
from observing H II regions older than a few million years. 

The continuous star formation cluster models often provide as good a fit to the observations 
as the single-age burst models. Although most published photoionization models are based on 
instantaneous burst models, this suggests that continuous models should also be considered as 
valid alternatives. Star formation extending over time (an age spread of a few Myr), as in the case 
of individual star clusters formed in subsequent bursts, but close enough to remain unresolved by 
the spectrograph slit, is a scenario that deserves further investigation. This is the case of 30 Dor, 
whose central ionizing cluster has an age of ~3 Myr; there is however evidence for previous episodes 
of star formation, as indicated by the presence of supergiants (Walborn 1991). These components of 
different age are spatially resolved, but they would not be if 30 Dor were at a much larger distance. 
Recent work (Garcia- Vargas etal. 1997, Gonzales-Delgado etal. 1997) indicates the coexistence of 
young (~3 Myr) and old (~8 Myr) components in individual starbursts, supporting the idea of an 
age spread, and suggesting star formation in multiple, subsequent bursts as a likely mechanism in 
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giant H II regions. 



5. TEMPERATURE OF THE IONIZING STARS 

Shields k Searle (1978) were the first to point out the importance of the [SIII]AA9069,9532 
lines to constrain the ionization structure of H II regions in the absence of temperature-sensitive 
lines, since sulfur is an effective coolant through these near-IR forbidden lines and the fine-structure 
lines at 19 and 33 fim. Mathis (1982, 1985) later introduced a method to determine the relative 
effective temperature scale, based on S + /S ++ and + /0. The use of ratios of ionic stages of the 
same element makes this method almost independent of nebular abundance. However, absolute 
values of T* are much more difficult to obtain, due to the dependence on the stellar atmospheres 
adopted for the nebular models. Besides, there is a saturation effect for very hot stars, so that the 
method becomes less useful at high temperatures (Garnett 1989). 

A 'radiation softness' parameter r\ was introduced by Vflchez k Pagel (1988) as a modification 
of Mathis' procedure, and defined as 

o + /o ++ (f{ . 

It is almost independent of ionization parameter, reddening, electron density and temperature. The 
observed quantity is given by 

_ [OII]AA3726,3729 / [OIII]AA4959, 5007 
V ' ~ [SII]AA6717,6731 / [SIII]AA9069, 9532. ^ ' 

As a measure of the hardness of the radiation field, this parameter provides us with an indicator 
of the %g of the exciting stars of H II regions. The large difference in the ionization potentials of + 
(35.1 eV) and S + (23.2 eV) is responsible for the sensitivity of 77 to the spectral energy distribution 
of the ionizing radiation. Spatially resolved observations of the Orion Nebula, NGC 604 in M33 
and 30 Doradus in the LMC confirm the reliability of this method (Vflchez k Pagel 1988, Mathis 
k Rosa 1991), although Ali etal. (1991) questioned it on the basis of possible effects of dielectronic 
recombination on the sulfur lines (the relative coefficient has not been computed yet). Applications 
of the 77 method to study the physical properties of extragalactic H II regions include the works of 
Diaz etal. (1987), Vflchez k Pagel (1988), Diaz etal. (1991) on the high-metallicity H II regions in 
M51, and Garnett k Kennicutt (1998) on a large sample of objects in M101. 

There is a general agreement about the existence of a metallicity dependence of the hardness 
of the radiation for H II region ionizing clusters. For example, in a sample of high-excitation H II 
regions Stasihska (1980) found that %jj increases statistically with decreasing abundance. The 
same conclusion has been reached by studies of the generally metal-poor H II galaxies (Campbell 
etal. 1986, Melnick 1992, Cerviho k Mas-Hesse 1994). Vflchez k Pagel (1988) confirmed this trend 
based on a compilation of H II region data from several authors. At any given abundance the 
spread in %g typically amounts to ~5,000 K, which can be due to several combined causes, such 
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as evolutionary effects, uncertainties in the models, and differences intrinsic to the ionizing popula- 
tion. The methods used in these works were based on the determination of electron temperatures, 
mainly via the [OIII]A4363 emission line, and are henceforth limited to low abundance (generally 
Z < 0.4 Z Q ). 

Different authors have instead debated over the cause of the T* gradient, whether it is a 
reflection of a change in the upper IMF (e.g. Campbell etal. 1986), or simply an effect of stellar 
evolution at different metallicities (McGaugh 1991). To illustrate the latter, we show in Fig. || the 
metallicity-2^ relation predicted by the Geneva stellar evolutionary models for different stellar 
masses, and for ages of zero and one Myr (Schaerer etal. 1993b, and references therein). For 
massive stars a decrease of the order of 5,000 K or more, up to ~10,000 K, exists in the metallicity 
range 0.05 Z Q - 2 Z . The effect is stronger at low abundance, where the gradient is steeper, which 
seems to be in agreement with H II region observations. 

In the following we examine our data in search for variations of T* with abundance, by means 
of the analysis of the rjf parameter. Since our sample covers the range in metallicity from about 
0.2 Zq to above solar, we are able to extend previous work in order to check if the T* trend continues 
at higher abundances. However, as mentioned above, the expected change in this metallicity range 
is rather limited, making it more difficult to detect. The fact that at high abundance the modelling 
of H II regions becomes problematic (due also to the lack of electron temperatures, forcing us 
to use an empirical calibration for the metallicity) complicates the interpretation of the data. In 
addition, the dependence of the results on the stellar atmospheres adopted, and the sensitivity to 
the ionization parameter make the precise T* scale, and in part its gradient, model dependent. 

5.1. gradient from single-star models 

Figure [?] shows the rjf vs R23 diagnostic diagram for our data and single-star models, and 
illustrates the main difficulties in deriving a T* gradient from rjf. The main complications are 
its dependence on R23 (at constant T and U) and the weaker sensitivity to T# in the high-T* 
range. Changes in the ionization parameter also have considerable effects on the results. We finally 
remind the reader about the rather mild T* dependence on metallicity which is expected when 
the abundance is larger than about half solar. The H II region parent galaxies are labeled in the 
rjf - R23 plot of Figure || where the sample has been divided into early-type (Sa-Sb; upper panel) 
and late-type (Sbc-Sm; lower panel) spirals. We also include spatially resolved observations of the 
30 Doradus nebula in the LMC by Kennicutt & French (Kennicutt et al, in preparation), which 
show that the dispersion in line ratios for a typical supergiant H II region is roughly of the same 
magnitude as the dispersion in the whole sample. When compared with the models of figure [?], 
this leads us to conclude that accurate T* determinations cannot be obtained through rjf. The 
strength of the method lies in the possibility of detecting large-scale trends of T* in a large sample 
of high- abundance H II regions. 
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In order to simplify the interpretation of the data, in Fig. ^a we have divided the H II region 
data into two ionization parameter classes (log U above and below —2.75), according to the mea- 
sured [SII] / [SIII] ratio and equation (Q), and superposed the models calculated for logU = —2.5 
and —3.0 at various stellar effective temperatures for the Kurucz atmospheres. The comparison of 
the models with the observed sequence suggests a metallicity dependence of T*, in the sense that 
the upper bound of T* shifts to lower temperatures at higher abundances, in a fashion similar to 
what Campbell (1988) noted for metal poor H II galaxies. Within our sample this effect is more 
easily seen in the high ionization parameter range (logfJ > —2.75). Quantifying it is however 
much more difficult, and, as the bottom part of Fig. ||a shows, the models do not agree with all 
the observational data. This discrepancy was seen by, among others, Skillman (1989) and Garnett 
(1989), who pointed out that the S + /S ++ ratio is underpredicted (or + /0 ++ overpredicted) by 
nebular photoionization models. It has been suggested that the origin of the disagreement lies in the 
sulphur atomic parameters (Garnett 1989) or in density inhomogeneities of the nebular gas (Rosa 
1993). Garnett (1989), after examining the effects of different stellar atmospheres on the model 
predictions for these two ionic ratios, concluded that only relative changes of T* can be inferred 
from 7]. Before reaching any conclusion on the existence of a T* gradient at high abundance we 
will also briefly consider the effects of changing the input stellar atmospheres, including the most 
recent, state-of-the-art hot star atmosphere models. 



5.2. Effects of input stellar atmospheres on the interpretation of the T* gradient 

The choice of stellar atmospheres, such as the widely used models by Mihalas (1972) and 
Kurucz (1979, 1992), is responsible for well-known differences in the predicted emission line ratios 
when computing H II region photoionization models (Skillman 1989, Evans 1991). These models 
have well-recognized inadequacies, namely the absence of line blanketing (Mihalas) or treatment of 
plane-parallel atmospheres in LTE (Kurucz). The inadequacy of these models for early- type stars 
derives from severe departures from LTE and the onset of radiation-driven stellar winds (Kudritzki 
& Hummer 1990, Kudritzki 1998). Enormous progress has been made in the last decade in stellar 
atmosphere modeling, in being able to treat theoretically the effects of departures from LTE, opacity 
from metal lines, stellar winds, and spherical extension on the computed hot star spectral energy 
distribution and ionizing flux (Gabler etal. 1989, Najarro etal. 1996, Schaerer etal. 1996a, b). The 
modifications with respect to previous calculations are profound. For instance, the predicted flatter 
ionizing flux distribution in the He I continuum (below 504 A) has already offered a solution for 
a long-lasting discrepancy between observations and nebular models regarding the intensity of the 
[Nelll] lines, severely underpredicted (up to a factor of 10 or more) when adopting the Kurucz 
atmospheres (Sellmaier etal. 1996, Stasihska & Schaerer 1997). Our understanding of the stellar 
populations responsible for the ionization of extragalactic giant H II regions is also affected, as 
shown by the possible explanation of high-excitation nebulae in terms of 'normal' massive O stars, 
instead of more exotic ionizing sources (Gabler etal. 1992). 



-16- 



Various predictions of the spectral energy distribution of hot stars are becoming available 
(Pauldrach etal. 1994, Sellmaier etal. 1996, Najarro etal. 1996, Schaerer & de Koter 1997, Paul- 
drach etal. 1998). While basically the same physical processes are taken into account, considerable 
differences exist in modeling techniques, for example in computing the occupation numbers for 
the different ions, so that different H II regions ionization structures would be predicted. While 
a detailed comparison is out of the scope of the present paper, we have computed nebular models 
based on stellar atmospheres different from Kurucz's, in order to verify that the general, qualitative 
conclusions remain unchanged. We have therefore computed nebular models using the most recent 
early-type star atmospheres available: the CoStar models (Schaerer & de Koter 1997), available for 
two different metallicities (Z and 0.2 Z ), and recent hydrodynamical non-LTE calculations of ex- 
panding atmospheres kindly made available to us by the Munich group (see Pauldrach etal. 1998). 
The latter, too, were calculated for Z = Zq, 0.2 Z . The 2 Z nebular models were computed 
with solar abundance stellar atmospheres. Additional models were computed with the Mihalas 
(1972) atmospheres, keeping the abundance of the stellar ionizing source fixed. All the remaining 
parameters were the same as described earlier for the Kurucz atmospheres. 

Without discussing the different theoretical approaches, for which we refer the reader to the 
above-mentioned papers, or the apparent merits and faults of the atmosphere models used in this 
comparison, we present in Figs. ^b-|^d the resulting models in the rjf - R23 diagram, with the same 
partition of the data into two ionization parameter bins as in Fig. |a. The fi gures demonstrate 
very clearly that the absolute T* scale greatly depends on the adopted input ionizing flux. It is 
particularly interesting that the two most recent non-LTE expanding atmosphere models give very 
different results in the lower %g regime; future work on H II region modeling can in principle be 
very useful in solving these discrepancies. Concentrating however on the points in common, we find 
that all models are consistent with higher T* values (> 50, 000 K) at the lowest abundances, while 
at high abundance T* levels off at a smaller value, between 35,000 K and 45,000 K, depending on 
the model atmospheres. We conclude that the data cannot be interpreted with a constant effective 
temperature at lower abundance: somewhere between Z and 0.2 Z (we cannot specify a precise 
value with the present models and observations) the average r)l— R23 relation steepens more than 
the constant temperature models. Constant T* cannot be ruled out for higher abundances. 



5.3. Trends at low and high abundance 

The line ratios studied earlier seem to conspire against a clear definition of a temperature 
dependence on metallicity. But, as seen at the beginning of this section, direct electron temperature 
determinations in nebular spectra below Z ~ 0.4 Z have led to the notion of a T^g gradient for 
metal poor objects. How can this be reconciled with the present work? In the diagnostic diagram 



of Figure 10 (top) the open stars mark the position of a few galactic and Magellanic Cloud H II 
regions with metallicity 0.2 < Z/Z < 0.5 from the work of Dennefeld & Stasihska (1983). These 
are rather low luminosity H II regions when compared with the objects in our sample, but they 
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might serve for the purpose of better understanding the diagnostic diagrams. The models shown 
for comparison were computed with the Munich group atmospheres at log U = —2.0, since these 
are high ionization parameter objects. No metallicity effect on T* can be derived from this plot. 
For comparison the bottom part of Figure |l0| shows the objects in our H II region sample for which 
an abundance based on electron temperature measurements is available in the literature. Here the 
models refer to logU = —2.5, in agreement with the average <logC/> = —2.6 for the data points. 
For these two sets of H II regions we compare in Fig. 11 the observed rjf - O/H relation with the 
one predicted by nebular models (Munich group atmospheres, logU = —2.5). The 30 Dor position 
is marked by the open star. The T* trend with metallicity becomes now clearer, and well agrees 
with previous work (see for example Fig. 6 in Stasihska 1980, or Fig. 4.3.2 in Melnick 1992). 

What can be said at the opposite (metal-rich) end of the abundance range? As mentioned 
earlier, the rjf data can be consistent with a constant T* above the solar metallicity. This statement 
can be checked by examining another important diagnostic at low temperatures, the intensity of the 
neutral helium line He I A5876 with respect to the H/3 emission line intensity. This line ratio provides 
a particularly sensitive diagnostic in the 35,000-40,000 K range, and tends to saturate at higher 
temperatures; unlike rjf, it is only weakly dependent on ionization parameter. The corresponding 
observations and nebular models are shown in Fig. 12, which suggests that H II regions with above- 
solar abundance define an upper bound in temperature at T* ~ 35, 000 K. Comparison of our models 
with the data strongly suggests that T* drops sharply below 40,000 K for Z > Z©. We regard this 
result to be quantitatively more accurate than the one obtained with rjf, given the uncertainties 
already described for the latter. Qualitatively, there is no contradiction between the two results. 

In conclusion, the single-star models support the idea of an important dependence of on 
metallicity over the whole range of abundance displayed in our sample. However, we must use some 
caution in the interpretation of this result, which is dependent on the various assumptions made 
in our nebular models, besides the uncertainties in the input ionizing fluxes. 



6. DISCUSSION: DOES THE IMF CHANGE WITH METALLICITY ? 

Investigations of possible variations of the IMF in star forming regions have often been stim- 
ulated by the Shields and Tinsley (1976) suggestion of a dependence of the upper mass limit on 
metallicity, M up oc Z~ a / 2 , with a ~ 1. This followed from theoretical considerations by Kahn 
(1974), who argued that the increase of the interstellar grain opacity with Z would, as an effect of 
the consequent increase in radiation pressure, prevent further accretion onto a forming star, and 
set an upper limit to the luminosity-to-mass ratio. We note, however, that Shields and Tinsley's 
theoretical relation between EW(H/3) and T u (temperature of the hottest stars present) was, con- 
sidering the measured gradient of EW(H/3) in M101, only marginally consistent with a dependence 
of M up on Z. Adopting a Salpeter (1955) slope (a = —2.35) for the IMF, instead of the steeper 
value used by Shields and Tinsley (a ~ —4), makes the inferred T u gradient inconsistent with a 
significant dependence of M. up on metallicity (Scalo 1986). 
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Direct counts of massive stars in resolved star clusters and OB associations in our Galaxy 
and the Magellanic Clouds do not show a variation of the IMF with the galactic environment or 
the metallicity, at least for metallicities equal to solar or below (Massey etal. 1995a, b, Hunter 
etal. 1997). For more distant star forming regions, studies of the nebular recombination lines {e.g. 
Garcia- Vargas etal. 1995, Stasihska & Leitherer 1996) and of the UV resonance lines (Robert, 
Leitherer & Heckman 1993) have also pointed to a 'universality' of the IMF. In this work we have 
compared emission-line diagnostics for a large H II region sample to theoretical predictions. We 
now conclude with some comments on our results, and their implications for constraints on the 
upper IMF variations. 

Our cluster photoionization models often are difficult to interpret unambiguously, due to the 
fact that different combinations of age, ionization parameter and IMF can reproduce the observed 
sequences in the diagnostic diagrams, as seen in Fig. |E[ Ambiguities are also introduced by the 
sensitivity of the inferred stellar temperatures to the choice of stellar atmosphere models. A further 



example is presented in Fig. 13 for the r/f parameter. The two IMFs with slope a = 2.35 and 
a = 3.30, but having the same upper mass limit 100 Mq, are virtually indistinguishable from a 
comparison of the observations with our models. The reason of this insensitivity of optical nebular 
spectra to the slope of the mass function is the relatively small fractional range in stellar mass 
among the ionizing stars. 

For M up = 30 Mq the predictions are somewhat different. The scatter of the data makes 
it difficult to say which IMF better represent the data, but we can exclude such a small upper 
mass limit for low-metallicity H II regions. To test the suggestion that the upper mass limit of 
the IMF is lower at higher metallicity we need to examine the abundance range Z > Zq. Even 
though in general the differences using the two upper mass limits are small, and probably within 
the dispersion of the data, there is one outstanding case, where the predictions are widely divergent: 
the He I A5876 emission line vs R23 diagram (Fig. |l4|). The low He I A5876/H/3 values measured at 
high abundance (~ 2 Zq) are in better agreement with the lower IMF mass cutoff, while the 100 Mq 
upper mass limit models still provide a good fit to the data at Z ~ Zq. This is consistent with the 
single-star results discussed in § || where an upper limit for T* around 35,000 K was established. 
This roughly corresponds to an 09V spectral type star, having a mass of 25 Mq (Vacca, Garmany 
& Shull 1996). An upper limit for T* of 40,000, which seems to be excluded, corresponds to a 
~ 35 Mq stellar mass. Correcting the He I emission lines by different, reasonable amounts of 
underlying absorption (equivalent width between zero and one) does not modify our conclusion 
substantially. The overall variation of T*, of the order of 15,000 K in the 0.2-2 Zq abundance 
range, inferred from our radiation hardness parameters, seems too large to be accounted for with 
a single upper mass limit for star formation, considering the current models of stellar evolution 
at different metallicity (see Fig. ||). If our understanding of stellar evolution is correct, our data 
and models therefore suggest an IMF biased toward a significantly lower upper mass limit for 
high-met allicity (between Zq and 2 Zq) H II regions. Whether this is an effect of a continuous 
metallicity gradient, as proposed by Shields and Tinsley (1976), or a more or less abrupt change in 
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the star formation properties at high abundance, cannot be judged from our work. Our conclusion 
is not necessarily in contrast with previous findings. Most authors have concentrated on rather low 
abundance H II regions, nevertheless the change with metallicity that we infer from their works 
is consistent with ours (see for example Melnick 1992). A constant upper IMF is still consistent 
with our data for Z < Z©. We must consider this result with some caution, since it is of course 
not only affected by the details of the nebular models (gas physical properties, photoionization 
code and atomic parameters), especially in the critical high-abundance regime, but also by the 
shape of the ionizing continuum (determined by the stellar evolutionary tracks and/or the model 
atmospheres). Uncertainties in any of these will propagate to the final result. The adoption of 
different stellar tracks, for example, can produce non-negligible effects on the computed emission- 
line intensities (Garcia- Vargas 1996, Stasihska & Leitherer 1996). The effects of the most recent 
atmosphere models, including departures from LTE and stellar winds, on the ionization structure 
of H II regions, need to be investigated in detail as well. 

In summary, our results imply that either: 1) The upper mass limit decreases with increasing 
metal abundance; or 2) Current stellar models underestimate the decrease in stellar temperature 
with increasing abundance, at a given mass. 

7. SUMMARY 

In this paper we have presented a new spectroscopic survey of 96 H II regions contained in 
spiral galaxies of different Hubble type. To interpret the measured nebular emission line ratios, two 
large sets of photoionization models have been calculated: single-star models based on the LTE 
Kurucz stellar atmospheres, and cluster models in which the ionizing continuum is taken from the 
evolutionary synthesis models of Leitherer & Heckman (1995). Additional single-star models were 
computed with more recent non-LTE stellar atmosphere models. 

We have mostly focused on two line diagnostics, the softness parameter rjf and the He I A5876/H/3 
line ratio, in order to investigate how the temperature of the ionizing stars depends on metallicity 
in our sample of high-abundance H II regions. From the single-star models we derive a rather 
important dependence of the upper bound of T* on metallicity, from 50,000 K or more at the 
low-abundance end (~ 0.2 Z©) to 35,000 K at above-solar abundances. With our current under- 
standing of stellar evolution at different metallicity, this suggests a smaller upper mass limit for 
star formation at high metallicity, as confirmed also by the cluster models. 

We have found that, interpreting H II regions as ionized by single-age populations of stars, only 
young (1 and 2 Myr) burst models can reproduce the observed emission line sequences. This implies 
either that the importance of the W-R phase is overestimated by the stellar evolution models, or 
that some physical mechanism disrupts the H II regions in a few million years. We have shown 
that continuous star formation models up to a few Myr of age can be well compared with the data, 
suggesting that in real H II regions populations of stars of different age can coexist. 
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The ionization parameter is found to vary approximately in the range —3.5 < log U < —2.5. 
Evidence for a slight trend of decreasing U with increasing abundance is found. Finally, an 
abundance-dependent N/O ratio is necessary to match the [Nil] /[Oil] model sequences to the 
observations. 

We thank G. Ferland for the existence of CLOUDY. We are indebted with A. Pauldrach 
and R. Kudritzki for making their stellar atmosphere models available to us. Comments from an 
anonymous referee helped us to improve the manuscript. 
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Table 1. List of abbreviations 



abbreviation meaning 

EW equivalent width 

IMF initial mass function 

a slope of a power-law IMF (dN/dm oc m~ a ) 

M-low IMF lower mass limit 

M up IMF upper mass limit 

Quo number of H ionizing photons/sec 

Qjj e o number of He ionizing photons/sec 

0+/Q++ 
V S+/S++ 

f [OII]AA3726,3729 / [QUI] AA4959, 5007 
^ [SII]AA6717,6731 / [SIII]AA9069,9532 

R 23 ([OII]A3727 + [OIII]AA4959,5007)/H/? 

SED spectral energy distribution 

U ionization parameter 
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rp i i o 

Table 2. 


The galaxy sample 






Galaxy 


Type a 




v a 


Source of 










(km s _1 ) 


Blue Spectra b 


NGC 224 


(M31) 


Sbl-II 


-21.61 


-10 


90" 


NGC 598 


(M33) 


Sc(s)II-III 


-19.07 


69 


90" 


NGC 628 


(M74) 


Sc(s)I 


-20.87 


861 


MRS 


NGC 1569 




SmIV 


-15.34 


144 


90" 


NGC 2403 




Sc(s)III 


-19.47 


299 


90", MRS 


NGC 2841 




Sb 


-20.65 


714 


MMT 


NGC 3031 


(M81) 


Sb(r)I-II 


-20.75 


124 


MMT 


NGC 5194 


(M51) 


Sbc(s)I-II 


-20.72 


541 


90" 


NGC 3310 




Sbc(r) 


-19.94 


1073 


90" 


NGC 3351 


(M95) 


SBb(r)II 


-19.78 


641 


MMT 


NGC 3368 


(M96) 


Sab(s)II 


-20.53 


758 


MMT 


NGC 3521 




Sbc(s)II 


-20.50 


627 


ZKH 


NGC 3623 


(M65) 


Sa(s)II 


-20.60 


675 


MMT 


NGC 4258 


(M106) 


Sb(s)II 


-21.17 


520 


MMT 


NGC 4736 


(M94) 


RSab(s) 


-19.93 


345 


MMT 


NGC 4861 




SBmlll 


-17.76 


836 


90" 


NGC 5236 


(M83) 


SBc(s)II 


-20.24 


275 


D80 


NGC 5701 




(PR)SBa 


-19.62 


1424 


MMT 


NGC 6384 




Sb(r)1.2 


-21.40 


1735 


MMT 


NGC 7331 




Sb(rs)I-II 


-21.72 


1114 


MMT 



a From Sandage & Tammann 1987; vo corrected to Ho = 75 km s 1 Mpc 1 

b 90": B&C + Steward Observatory 90-inch telescope; MMT: Red Channel + MMT; MRS: 
McCall etal. 1985; ZKH: Zaritsky etal. 1994; D80: Dufour etal. 1980 
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TABLE 3 
Reddening-corrected line fluxes 



No. 


Offs 




[Oil] 


[OIII] 


Hel 


[Nil] 


[SII] 


[SHI] 


Name b 




RA 


DKC 


3727 


49594-5007 


5876 


65484-6583 


67164-6731 


90694-9532 




M31 


1 


-157 


-3244 


209 ± 11 


347 ± 14 


11 ± 1 


22 ± 1 


11 ± 1 


35 ± 5 


P7 


2 


-840 


-1498 


160 ± 8 


40 ± 2 


13 ± 3 


127 ± 5 


47 ± 2 


86 ± 14 


P258 


3 


-1165 


-717 


184 ± 14 


39 ± 11 


10 ± 5 


88 ± 8 


66 ± 6 


50 ± 8 


P286 


4 


-838 


-663 


206 ± 11 


30 ± 1 


12 ± 1 


174 ± 7 


81 ± 3 


71 ± 11 


P309 


5 


561 


-377 


223 ± 11 


87 ± 3 


16 ± 3 


137 ± 6 


66 ± 2 


70 ± 10 


P420 


6 


867 


1972 


330 ± 18 


86 ± 4 


35 ± 2 


89 ± 4 


43 ± 1 


49 ± 7 


P729 


M33 


1 


-824 


-24 


156 ± 8 


641 ± 26 


13 ± 1 


21 ± 1 


25 ± 1 


95 ± 15 


NGC 588 


2 


553 


420 


224 ± 11 


235 ± 9 


13 ± 1 


32 ± 2 


39 ± 1 


95 ± 15 


NGC 604 


3 


-220 


113 


237 ± 12 


169 ± 7 


15 ± 1 


55 ± 2 


24 ± 1 


124 ± 19 


NGC 595 


M51 


1 


34 


12 


38 ± 10 


7 ± 8 


± 3 


74 ± 6 


24 ± 6 


± 


CCM 37a 


2 


15 


30 


± 10 


± 10 


± 5 


70 ± 8 


21 ± 9 


± 


CCM 6A 


3 


3 


57 


27 ± 3 


± 1 


9 ± 1 


105 ± 4 


40 ± 2 


20 ± 4 


CCM 107 


4 


70 


138 


106 ± 6 


18 ± 3 


10 ± 2 


162 ± 6 


47 ± 3 


53 ± 7 


CCM 10 


5 


-83 


-108 


70 ± 6 


± 4 


± 2 


121 ± 7 


49 ± 7 


46 ± 9 


CCM 71 


6 


-88 


-82 


57 ± 3 


8 ± 1 


9 ± 1 


150 ± 6 


49 ± 2 


60 ± 9 


CCM 72 


7 


90 


-135 


88 ± 5 


22 ± 3 


9 ± 1 


153 ± 6 


42 ± 3 


65 ± 10 


CCM 55 


8 


104 


-116 


104 ± 6 


24 ± 3 


10 ± 2 


162 ± 6 


46 ± 3 


77 ± 12 


CCM 53 


9 


-66 


50 


83 ± 17 


35 ± 21 


9 ± 8 


172 ± 13 


45 ± 11 


98 ± 12 


CCM 84 


10 


-60 


119 


157 ± 16 


42 ± 13 


± 6 


144 ± 13 


47 ± 10 


± 


CCM 91 


11 


-134 


-188 


167 ± 10 


70 ± 5 


10 ± 3 


161 ± 8 


70 ± 4 


97 ± 16 


CCM 71A 


12 


-59 


-125 


91 ± 9 


7 ± 7 


± 3 


168 ± 8 


54 ± 6 


± 


CCM 68 


13 


-23 


-122 


48 ± 7 


± 4 


± 2 


120 ± 7 


45 ± 7 


± 


CCM 64 


M81 


1 


107 


-291 


249 ± 27 


101 ± 14 


13 ± 6 


115 ± 9 


36 ± 7 


81 ± 13 


GSl 


2 


40 


-250 


262 ± 38 


121 ± 16 


13 ± 7 


111 ± 9 


41 ± 8 


63 ± 13 


GS2 


3 


-290 


89 


339 ± 27 


252 ± 19 


13 ± 5 


93 ± 8 


30 ± 6 


115 ± 17 


GS4 


4 


-273 


-19 


340 ± 27 


106 ± 14 


12 ± 7 


135 ± 11 


71 ± 10 


81 ± 14 


GS7 


5 


-196 


292 


222 ± 22 


89 ± 13 


13 ± 6 


159 ± 10 


80 ± 8 


77 ± 12 


GS9 


6 


-92 


300 


191 ± 22 


159 ± 17 


13 ± 7 


122 ± 10 


49 ± 8 


67 ± 11 


GSll 


7 


-47 


260 


334 ± 39 


71 ± 26 


12 ± 9 


126 ± 12 


51 ± 11 


83 ± 13 


GS12 


8 


146 


7 


139 ± 33 


48 ± 31 


12 ± 16 


141 ± 21 


41 ± 20 


68 ± 16 


GS13 


9 


237 


-845 


209 ± 36 


902 ± 73 


14 ± 10 


48 ± 13 


34 ± 13 


54 ± 12 


Miinchl 


10 


-142 


526 


156 ± 31 


415 ± 23 


15 ± 5 


76 ± 6 


36 ± 6 


76 ± 11 


Munchl8 


NGC 628 


1 


49 


52 


196 ± 16 


48 ± 10 


14 ± 5 


128 ± 8 


38 ± 6 


90 ± 14 


H292 


2 


-74 


-22 


152 ± 15 


23 ± 13 


10 ± 7 


114 ± 10 


56 ± 8 


58 ± 9 


H451 


3 


-60 


-107 


245 ± 18 


37 ± 6 


11 ± 3 


129 ± 7 


64 ± 4 


81 ± 13 


H572 


4 


42 


-116 


255 ± 22 


58 ± 18 


12 ± 10 


104 ± 13 


54 ± 12 


72 ± 12 


H598 


5 


-42 


-154 


235 ± 15 


189 ± 14 


14 ± 6 


53 ± 8 


30 ± 8 


90 ± 14 


H627 


6 


-186 


86 


418 ± 23 


111 ± 8 


12 ± 3 


81 ± 5 


54 ± 4 


84 ± 13 


H154-155 


7 


-89 


192 


325 ± 17 


157 ± 8 


12 ± 2 


72 ± 4 


39 ± 3 


71 ± 10 


H13 


NGC 1569 


1 


-7 


3 


116 ± 6 


922 ± 37 


20 ± 1 


6 ± 


12 ± 


50 ± 7 


HK45 


2 


25 


-12 


100 ± 6 


895 ± 36 


19 ± 1 


10 ± 1 


18 ± 1 


55 ± 5 


HK23 


NGC 2403 


1 


-96 


35 


300 ± 15 


257 ± 10 


20 ± 2 


80 ± 4 


75 ± 3 


118 ± 19 




2 


-45 


55 


210 ± 11 


148 ± 6 


14 ± 1 


74 ± 4 


47 ± 3 


111 ± 17 




3 


-186 


45 


241 ± 12 


236 ± 9 


14 ± 1 


44 ± 2 


39 ± 2 


98 ± 15 




4 


10 


32 


228 ± 12 


184 ± 7 


14 ± 1 


51 ± 2 


49 ± 2 


110 ± 16 




5 


-133 


-146 


258 ± 13 


460 ± 19 


13 ± 2 


30 ± 2 


41 ± 3 


111 ± 16 




6 


63 


-49 


234 ± 20 


35 ± 19 


9 ± 11 


116 ± 16 


86 ± 16 


68 ± 13 




7 


165 


136 


207 ± 12 


461 ± 24 


12 ± 4 


12 ± 6 


12 ± 6 


63 ± 9 




8 


45 


69 


365 ± 41 


117 ± 37 


14 ± 20 


50 ± 24 


28 ± 24 


88 ± 18 




9 


-494 


137 


406 ± 25 


235 ± 57 


10 ± 10 


24 ± 13 


39 ± 13 


53 ± 11 




NGC 2841 


1 


29 


-116 


134 ± 14 


22 ± 13 


16 ± 7 


154 ± 11 


61 ± 10 


± 


HK27 


2 


-42 


86 


64 ± 13 


15 ± 11 


± 5 


113 ± 9 


36 ± 8 


± 


HK51 


3 


-42 


77 


43 ± 7 


5 ± 7 


± 3 


127 ± 10 


61 ± 9 


± 


HK52 


4 





71 


87 ± 7 


9 ± 4 


± 2 


145 ± 9 


40 ± 7 


54 ± 10 


HK62 


NGC 3310 


1 


-9 


-6 


267 ± 14 


380 ± 16 


12 ± 1 


53 ± 2 


41 ± 1 


74 ± 12 


HK53 


2 


6 


1 


408 ± 21 


204 ± 9 


9 ± 2 


47 ± 4 


34 ± 3 


73 ± 12 


HK35 


3 


1 


5 


269 ± 14 


314 ± 13 


13 ± 1 


61 ± 3 


36 ± 1 


97 ± 16 


HK85 



-28- 



TABLE 3 — Continued 



No. 


Offs 




[Oil] 


[OIII] 


Hel 


[Nil] 


[SII] 


[SHI] 


Name b 




RA 


DKU 


3727 


49594-5007 


5876 


65484-6583 


67164-6731 


90694-9532 






NGC 3351 


1 


14 


67 


32 ± 4 


10 ± 4 


± 2 


77 ± 6 


30 ± 8 


± 


HK45 


2 


-58 


-135 


104 ± 6 


23 ± 3 


9 ± 2 


160 ± 7 


36 ± 3 


82 ± 14 


HK80 


3 


-88 


-26 


56 ± 4 


11 ± 3 


9 ± 2 


118 ± 5 


33 ± 3 


52 ± 10 


HK101 


4 


-105 


-86 


137 ± 7 


32 ± 2 


11 ± 2 


179 ± 7 


45 ± 3 


78 ± 13 


HK109 


5 


29 


-192 


228 ± 13 


93 ± 8 


9 ± 4 


142 ± 8 


59 ± 5 


75 ± 11 




NGC 3368 


1 


54 


-24 


81 ± 5 


18 ± 3 


8 ± 3 


151 ± 8 


83 ± 4 


100 ± 13 




2 


-52 


25 


131 ± 8 


24 ± 4 


± 1 


149 ± 9 


56 ± 8 


141 ± 25 














NGC 3521 










1 


13 


97 


189 ± 21 


110 ± 44 


14 ± 1 


109 ± 4 


47 ± 2 


82 ± 13 




2 


-5 


80 


119 ± 15 


20 ± 10 


14 ± 1 


150 ± 6 


49 ± 2 


60 ± 9 




3 


-48 


47 


220 ± 21 


37 ± 12 


13 ± 1 


123 ± 5 


57 ± 2 


51 ± 8 




NGC 3623 


1 


-3 


-86 


123 ± 6 


17 ± 5 


10 ± 3 


109 ± 6 


46 ± 6 


32 ± 7 




NGC 4258 


1 


-318 


183 


364 ± 19 


100 ± 6 


± 1 


97 ± 9 


47 ± 9 


± 




2 


-152 


-37 


333 ± 17 


209 ± 9 


11 ± 3 


91 ± 6 


87 ± 6 


71 ± 17 




3 


-32 


141 


236 ± 12 


92 ± 4 


± 1 


126 ± 7 


82 ± 6 


86 ± 16 




4 


-25 


130 


219 ± 11 


128 ± 5 


10 ± 2 


104 ± 4 


47 ± 3 


103 ± 16 




5 


-19 


121 


248 ± 13 


95 ± 4 


18 ± 2 


118 ± 5 


66 ± 4 


112 ± 17 




NGC 4736 


1 


37 


9 


206 ± 11 


117 ± 5 


12 ± 1 


134 ± 5 


48 ± 2 


103 ± 16 


HK8 


2 


37 


3 


199 ± 10 


67 ± 3 


13 ± 2 


175 ± 7 


56 ± 2 


121 ± 19 


HK9 


3 


-29 


14 


174 ± 9 


45 ± 4 


9 ± 2 


155 ± 7 


54 ± 3 


57 ± 8 


HK45 


4 


-38 


9 


173 ± 9 


34 ± 2 


9 ± 1 


174 ± 7 


63 ± 3 


67 ± 10 


HK52 


5 


-39 


20 


203 ± 10 


90 ± 4 


12 ± 1 


156 ± 6 


29 ± 1 


97 ± 14 


HK53 


NGC 4861 


1 








115 ± 6 


799 ± 33 


11 ± 1 


5 ± 1 


17 ± 1 


59 ± 9 




2 


6 


6 


192 ± 10 


653 ± 27 


12 ± 2 


9 ± 3 


26 ± 4 


68 ± 11 




NGC 5236 


1 


53 


58 


46 ± 12 


4 ± 13 


8 ± 7 


117 ± 11 


43 ± 8 


42 ± 7 


D2 


2 


-103 


11 


95 ± 16 


5 ± 13 


8 ± 7 


126 ± 10 


51 ± 8 


26 ± 7 


D3 


3 


-224 


-27 


145 ± 18 


45 ± 13 


13 ± 7 


135 ± 10 


40 ± 8 


76 ± 12 


D5 


4 


-262 


36 


126 ± 16 


21 ± 13 


12 ± 6 


153 ± 10 


54 ± 7 


58 ± 9 


D6 


NGC 5701 


1 


-7 


-87 


209 ± 11 


186 ± 9 


15 ± 2 


117 ± 5 


45 ± 5 


103 ± 15 




2 


79 


-60 


216 ± 12 


115 ± 6 


12 ± 4 


137 ± 10 


61 ± 8 


84 ± 16 




NGC 6384 


1 


-2 


-83 


181 ± 11 


50 ± 5 


± 2 


148 ± 8 


78 ± 6 


98 ± 14 




2 


51 


37 


81 ± 10 


5 ± 7 


± 4 


120 ± 27 


31 ± 13 


± 




3 


53 


45 


115 ± 10 


18 ± 8 


± 4 


145 ± 12 


110 ± 6 


67 ± 11 




4 


-32 


76 


266 ± 16 


71 ± 8 


17 ± 5 


122 ± 8 


72 ± 10 


± 




5 


-41 


70 


256 ± 15 


92 ± 8 


± 3 


95 ± 16 


50 ± 14 


± 




6 


130 


110 


285 ± 17 


225 ± 12 


10 ± 4 


72 ± 8 


102 ± 13 


84 ± 14 




NGC 7331 


1 


-26 


-52 


225 ± 21 


43 ± 13 


± 5 


129 ± 12 


52 ± 10 


± 


HK95 


2 


-31 


-27 


232 ± 16 


38 ± 8 


± 3 


149 ± 9 


105 ± 9 


± 


HK98 


3 


-32 


16 


200 ± 21 


42 ± 17 


± 5 


143 ± 10 


69 ± 11 


± 


HK101 


4 


27 


80 


250 ± 16 


126 ± 11 


± 3 


130 ± 8 


75 ± 7 


62 ± 12 


HK40 



NOTE. — Fluxes given in units of H/3 = 100 

a Relative to galaxy nucleus, positive in E and N directions, in arcseconds 

b Meaning of catalog names: HK: Hodge & Kennicutt 1983; H: Hodge 1976; P: Pellet etal. 1978; 
CCM: Carranza etal. 1969; GS: Garnett & Shields 1987; D: Dufour etal. 1980 
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Table 4. Adopted solar abundances 



Element log abundance a 

He -1.00 

C -3.44 

N -3.95 

O -3.07 

Ne -3.91 

Mg -5.42 

Al -6.53 

Si -4.75 

S -4.79 

Ar -5.44 

Ca -6.64 

Fe -5.33 

Ni -6.75 

Na -6.67 



relative to H 



Table 5. Qhc /Qh° from model stellar atmospheres 



%(K) 2 Z Z 0.25 Z 0.1 Z 



60,000 0.35 0.38 0.45 0.48 

50,000 0.26 0.28 0.32 0.34 

45,000 0.20 0.22 0.25 0.26 

40,000 0.10 0.10 0.10 0.11 

35,000 0.01 0.01 0.01 0.01 
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Table 6. Qh^'/Qh" from cluster models 



age (Myr) 2 Z Z 0.25 Z Q 0.1 Z, 



M u „ = 100 a = 2.35 



1 




0.19 


0.22 


0.30 


0.35 


2 




0.11 


0.15 


0.24 


0.29 


3 . . 




0.27 


0.22 


0.13 


0.16 


4 




0.29 


0.21 


0.18 


0.10 


5 . . 




0.29 


0.31 


0.14 


0.05 


6 . . 




0.24 


0.25 


0.01 


0.03 


i\/r 


on 


a = z.oo 








1 




0.07 


0.08 


0.12 


0.15 


2 




0.05 


0.07 


0.11 


0.14 


3 




0.03 


0.05 


0.09 


0.13 


4 




0.01 


0.03 


0.06 


0.11 


5 




0.18 


0.01 


0.03 


0.06 


6 




0.24 


0.20 


0.01 


0.03 




= 100 


a = 3.30 








1 




0.17 


0.19 


0.26 


0.30 


2 




0.10 


0.13 


0.21 


0.26 


3 




0.18 


0.15 


0.12 


0.16 


4 




0.23 


0.15 


0.13 


0.10 


5 




0.25 


0.25 


0.10 


0.05 


6 




0.20 


0.21 


0.01 


0.03 



-31 - 



0.5 



CM 



o 



o - 



C7> 

q -0.5 



~ i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r 



~i — i — i — I — i — i — r 



o o 
o o 

•o°o 
c • • 

<? o • °° 
o • 



o o »fcR) •# 

q o °» * ° °. r 



o • • 



I I I 



-I I I I I L 



-J I I I I 1_ 



-1.5 -1 -0.5 

Log ([Olll] 5007/H/?) 



0.5 



\ 

o 
o 
in 



1 h 

0.5 
h 



,0 -0.5 \r 



-1.5 - 



-2 



T 



• ra. 



*1 _n *» O 



OO 



O • CD 



2?- 



_i i i i_ 



1.5 -1 
Log ([Nil] 6583/Ha) 



-0.5 



Fig. 1.— Excitation sequences in terms of (a) [OII]/H/3 vs [OIII]/H/3 and (b) [OIII]/H/3 vs [Nil] /Ha. 
Filled circles indicate observations from this work, open circles are data from McCall etal. (1985). 
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1.5 2 2.5 

EW(H/8) (A) 

Fig. 2. — Comparison with the models of Stasihska & Leitherer (1996). The latter are shown for an 
ionizing cluster mass M = 10 3 and M = 10 6 M Q , with dotted and dashed lines, respectively. Our 
models for log U = —2 and log U = —3 (calculated for n = 10 cm~ 3 ) are shown by the continuous 
lines connecting triangle and square symbols, respectively. All models are for solar composition. 
(Top) [OIII]A5007/H/3 vs EW(H/3). (Bottom) [OII]A3727/H/3 vs EW(H/3). 
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Fig. 3. — Single-star model results for [SII] / [SIII] vs Log R23, for logU = —2 (top), logt/ = —3 
(middle) and log?7 = —4 (bottom). The lines connect points at three metallicity values (0.25 Z Q , 
Z and 2 Z Q ). The stellar %g is coded as shown in the top panel. The dots are the observed values. 
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Fig. 4. — The [Nil] /[Oil] vs R23 diagnostic diagram from models without (left) and with (right) a 
metallicity-dependent N/O ratio. Symbols as in Figure ||. 
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Fig. 5. — Cluster model results for [0II]/[0III] vs Log R23. The two top panels refer to ages 1 
and 2 Myr, the two lower panels to ages 3 and 4 Myr. Each of these panels gives results for the 
Salpeter (top), M up = 30 M (middle), and a = 3.30 (bottom) IMFs. The lines connect points at 
three metallicity values (0.25 Z , Z Q and 2 Z ). The ionization parameter is coded as shown in 
the first panel. The dots are the observed values. 
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Fig. 6. — The metallicity-2^ relation predicted by the Geneva stellar evolutionary models. The 
dotted lines refer to ZAMS models, the continuous line to 1 Myr models. The stellar mass is 
indicated on the left. 



-37- 




Fig. 7. — Single-star (Kurucz atmospheres) model results for r\i vs Log i?23- 
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Fig. 8. — Identification of the H II region parent galaxies in the rjf - R23 diagram. Early-type (Sa- 
Sb) and late-type (Sbc-Sm) spirals are shown in the upper and lower half of the figure, respectively. 
Spatially resolved observations of 30 Dor are included for comparison. 
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Fig. 9. — The nebular models based on the four theoretical stellar atmospheres discussed in the text 
are compared in these rjf vs Log R23 diagrams. For each of the four sets of photoionization models 
the data points shown for comparison have been divided according to the inferred value of log U, 
according to the [SII]/[SIII] ratio and the models based on the Kurucz atmospheres: log U > —2.75 
in the upper half, and log U < —2.75 in the lower half of each of the four panels. The corresponding 
nebular models are calculated for logU = — 2.5 and logU = —3.0, respectively. 
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Fig. 10. — (Top) The open star symbols mark the position in the 77/-.R23 diagram of a few galactic 
and Magellanic Cloud H II regions having 0.2 < Z/Z < 0.5 (from Dennefeld & Stasihska 1983). 
Models calculated for log U = —2.0. (Bottom) The objects in our sample that have O/H determined 
from electron temperature measurements are marked by the larger symbols in this diagram. Models 
calculated for log U = —2.5. 
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Fig. 11. — The T* dependence on abundance as determined using rff. The data included are 
galactic and Magellanic Cloud H II regions from Dennefeld Sz Stasihska (1983; triangles), H II 
regions from our extragalactic sample having O/H determinations (solid dots), and the average of 
spatially resolved observations of 30 Dor (Kennicutt etal. , in preparation; open star). The models 
superposed are based on the Munich group atmosphere models. 
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Fig. 12. — Observations and single-star (Munich group atmospheres) nebular models for the 
He I A5876/H/3 diagnostic. 
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Fig. 13. — Cluster model results for rji vs Log i?23- See Figure ||| for symbol explanations. 
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Fig. 14. — Observations and cluster nebular models calculated for an age of one (left) and two 
(right) Myr for the He I A5876/H/3 diagnostic. The top panels show models calculated with an 
upper mass limit M up = 100 M , the bottom panels for M up = 30 M . Symbols as in Figure ||| 



